The Effects of dj;2_y2-dxy Mixing on Vortex Structures and 

Magnetization 

Mahn-Soo Choi and Sung-Ik Lee 

Department of Physics, Pohang University of Science and Technology, Pohang 790-784 South 

Korea 
(February 1, 2008) 

Abstract 

The structure of an isolated single vortex and the vortex lattice, and the 

magnetization in a d-wave superconductor are investigated within a phe- 

nomenological Ginzburg-Landau (GL) model including the mixture of the 

dr^2_y2-wscve and daj^-wave symmetry. The isolated single vortex structure in 

a week magnetic field is studied both numerically and asymptotically. Near 

the upper critical field Hc2 , the vortex lattice structure and the magnetization 

are calculated analytically. 

PACS numbers: 74.20De, 74.60-w, 74.60.Ec 
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I. INTRODUCTION 



Since the discovery of the high temperature superconductors (HTSC's), a great number 
of experimental and theoretical investigations have been carried out to identify the symmetry 
of the superconducting pairing state. Although there still remain controversies, it is believed 
that the (ia;2_j^2-wave symmetry is most probable in HTSC's; e.g. a series of strong evidences 
have been provided by the phase sensitive experimentsH. 

In recent years, the structure of vortices in a d^2_y2 superconductor was of great interest. 
It was expected that the structure of d-wave vortices might be very different from that of 
s-wave vortices. And connected with the vortex structure, the question of order parameter 
symmetry admixture arose. 

From a standpoint of the spontaneous symmetry breaking, in a bulk with a perfect crystal 
symmetry, the order parameter transforms according to an irreducible representation! of 
the crystallographic point grou p&i. Thus in a crystal of tetragonal (-D4) symmetry, only 
one symmetry state (e.g. d^2_y2-w8ive) is allowed. However near interfaces, surfaces or 
impurities, the crystalline symmetry is not perfect and the dx2_y2-wave order parameter 
fluctuates spatially and hence induces components of other symmetry. This is also the case 
in the presence of vortices, to which we will confine ourselves in this work. 

Soininen et a/.i, starting from a simple microscopic model Hamiltonian and its 
Bogoliubov-de Gennes equation, found a substantial admixture of an induced s-wave and 
the dominant (ia;2_y2-wave. Based on this work, Berlinsky et a/.i and Franz et a/.i investi- 
gated the structure of rf-wave vortices using the Ginzburg-Landau(GL) theory. Ren et a/.i 
microscopically derived the GL equation from the Gorkov equation of a continuum mean 
field model, and it was studied in detail by Xu et al. for the (i-wave vortex structures. 
Joyntlli also included the d-s mixing term in their phenomenological GL model, from the 
symmetry argument. Very recently, Maki and Beal-MonodEl and Won and Maki0 also in- 
corporated d-s mixture. They used an interaction potential with the s channel repulsive 
Coulomb interaction in their weak-coupling model. In most of the works above, they ob- 
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served that the s-wave component induced near the vortex core causes the d-wave component 
to have fourfold anisotropy,lll0'0 and that the vortices arrange in an obhque lattice instead 
of triangular one.3li0'0 

Recently, the possibility of admixture of (ia;2_j^2-wave and d^iy-wave was also suggested. 
Koyama and Tachikilli took into account the internal orbital motion of the pairing elec- 
trons in their linearized Gorkov-type gap equation near Tc. They found that the dominant 
(ia,2_y2-wave symmetry was mixed with d^y-wave symmetry, instead of the s-wave. The cou- 
pling of the two (i-wave components led to a significant paramagnetic effect and to strong 
enhancement of the upper critical field at lower temperatures, which remarkably fits the 
experimental results on over-doped cuprate superconductors Bi2Sr2CuOy and Tl2Ba2Cu06. 
And Ichioka et a/.0 considered both drc2_y2-s mixing and d^2_y2-dxy mixing in the frame work 
of the classical Eilenberger equations. They studied the isolated single vortex structures and 
found that the amplitude of the (i^^y-wave component has the shape of an octofoil. 

As described above, for the case of d-s mixture, the single vortex and/or vortex lattice 
structure have been extensively studied0i'0, but not for the case of dx2_y2-dxy mixture. In 
this paper, we present analogous studies based on the GL theory which includes the mixing 
of (ia.2_y2-wave and d^'y-wave, assuming the dx2-y2-mme symmetry of the superconducting 
ground state. 



II. GL THEORY FOR A i:»-WAVE SUPERCONDUCTOR 



We start with the GL free energy functional proposed by Koyama and Tachiki: 



G{H) = J (dx) 
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where the order parameter (^I/-) corresponds to the (ia.2_y2-wave (drcy-wave) symmetry. 
As we previously assumed, in Meissner state, the only non-vanishing component is ip^ and 
a;-|-(T) < at T < Tc. Other coefficients are assumed to be positive and independent of 
temperature. 

In passing, we make a couple of remarks. In Eq. (|I|) we included the coupling term 
~ (\I'^^\E'^ + h.c), which was omitted in the model originally proposed by Koyama and 
Tachikilili. This term is necessary in order to get the most general (up to the fourth order) 
free energy functional from the symmetry consideration.! The mixing term ~ (\I'!^^_ + 
h.c.)Bz in the quadratic order give rise to the paramagnetic current.S As we will see below, 
it dominates other mixing terms and significantly affects the mixed state properties of (i-wave 
superconductors. 

In mean field approximation, neglecting thermal fluctuation effects, physical properties 
are described by the corresponding GL equations. It is convenient to introduce the dimen- 
sionless quantities by adopting the fundamental length scale to be the GL coherence length 
^+ of (See Table. |). In this case, the GL equations and Maxwell equations are written 



as 



nV+ - ^+ + |^+rv^+ + x+l^-rv^+ + C+r+i^t - y+B^i^- = o, (2) 
n2^_ + + T|V'_|'^- + x-\i^+?i^~ + C-^-V'+ - i^-B,ij+ = o, (3) 
= -l(^^n^^ + h.c.) - ^fi-\iij*_Uiij_ + h.c) 

+ ^z/+ V X (^;V^_ + h.c.)z, (4) 

where V x V x A = J and 11 = (— iV + A) and other parameters are defined in the Table 
|. In Eq. (^), the last term is the paramagnetic current contribution. 



III. ISOLATED SINGLE VORTEX STRUCTURE 

Consider an isolated single vortex near the lower critical field -ffci^^- For the problem 
of isolated single vortex, it would be convenient to decompose the order parameter into the 



form ip± = f± e'^'"'^^ . Since the rf^jy-wave component is induced through the direct couphng 
to the magnetic field, in the low field near Hd its amplitude is expected to be very small 
compared with dx^_y2-waMe component. Thus just for intuitive understanding, for a moment 
we neglect the effect of the coupling in the fourth order {x± — 0, C± — 0). The effect will 
be considered below. 

Let's first look at the phase distributions associated with the vortex. It is obvious from 
the rotational symmetry of the GL equations that the phases should have ip{r, 9)± = —9 up 
to an additive constant^, which we take to be zero; the two (i-wave components have the 
same winding. Thus, the differential equations are rewritten, in terms of f± only, as 



nV+ -/+ + /!- y+Bf_ = 



(5) 



nV- + H/_ + T/! - v.BU = 



(6) 



= -{fi + i^-'ft) (-1/r + A)0 + Z/+ V X (/+/_)£, 



(7) 



where = -|i - If + (-1 + A)\ A = A{r)e, B = B{r)z, and J = J{r)e. 

Second, near the vortex center (r — > 0), B{r) and f±{r) have the asymptotic behavior of 
the following form: 



where0 



B{r) = Bo + B^r^ + 0\r^] 
/+(r) = Cor(l + C2r2 + 0[r4; 



^0 = 
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(8) 
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(10) 



(11) 
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Note that since Bq ~ in the extreme type-II superconductors (k ^ 1) the asymptotic 

behaviors of B{r) and /+(r) deviates very httle from those of the conventional s-wave GL 
theory {C2 ~ -(1 + Bo)/S, B2 ~ -CI/2k^). 

Next, consider the region far from the vortex center (r ^ k). Assuming the extreme 
type-II superconductor (k ^ 1), we can treat /+ = 1 and neglect |V^/|_ ~ /-/ft^ <^ /_• 
Then by taking the curl of Eq. we get the usual London equation: 

K^V X V X B + B = 27r52(r), (14) 

and from Eq. (^) 

/+(r) ~ ^fi(r) ~ yK7^e-"/^ (r > /t). (15) 

Thus far outside the core, only the pure (ia;2_j^2-wave component remains. 

Now we discuss the effect of the mixing terms in the fourth order, i.e. with finite x± and 
(±. In general, due to the coupling of the form {ip^tp'^ + h.c.) the GL equations do not have 
spherical symmetry. However, within the assumption of small d^y-wave component, we can 
apply a partial wave expansion for 

^_ = E/^"^We^"'- (16) 

n 

It is straightforward to show that the only non- vanishing component is of = 1, and that 
the asymptotic results given above have no change. 

We also provide the numerical results for the distribution of the order parameter and 
the magnetic induction in Fig. |l|. As shown in the Fig. |l] (a), the magnitude f-{r) of the 
dxy-"w&ve component increases as the d^^-yi-dxy coupling strength z/+ (i/_) increases and as 
the temperature T decreases. However, /_(r) is so small that its effect on /+(r) and B{r) 
is negligible; the isolated single vortex structure of this model is very similar to that of the 
conventional s-wave GL theory. 
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IV. VORTEX LATTICE AND MAGNETIZATION NEAR THE UPPER 

CRITICAL FIELD 



In the vicinity of the upper critical field Hc2, the magnitudes of the order parameters 
are small and thus the non-linear terms in the GL equations Eqs. (0J^) are negligible. The 
magnetic field are also assumed to be constant all over the space since the inter-vortex 
spacing is much less than the penetration depth A+. Then the GL equations are reduced to 
a linearized form: 

(n^ - 1)^+ - z/+B,z^. = 0, (17) 
(tf + - U-B,ij+ = 0. (18) 

Each of the solutions ipL± is just a linear combination of the wave functions in the lowest 
Landau level, which is infinitely degenerate, and determined so as to minimize the GL free 
energy. 

We minimize the free energy by generalizing the Abrikosov's procedures^!'!!. In these 
procedures, the two components of the order parameter satisfy 

^ = LE^^^ (19) 
and that the current contribution is given by 

^^^Jl = X + - ^+{rL+'^L- + h.C.)]z. (20) 

Then the magnetization!! and the GL free energy are given by 



(22) 



(2«:2^ - l)f3A_ 

where Pa = (/+)/( /+)^ is a still undetermined parameter, and nls{T) = h? G{T) / F'^{T) 
with 



R{T) 



1_ H,2iT) - 1 



F{T) = 1 - 2u+R + fi-^R'^ 



G{T) = l + 2{x+ + C+)R' + fi-'rR\ 



The Abrikosov parameter concerns the vortex lattice structure and is determined by 
minimizing the GL free energy. Since has the same form as the s-wave Abrikosov 
solution, it is obvious that the free energy is minimized for the triangular vortex lattice 
(Pa = 1 .16)11. The magnetization in Eq. ( pT]) shows the usual form of the conventional 
s-wave GL theoryilil, but the slope AirdM/dH depends on the temperature. 



There is no direct observation of the single vortex structure in a HTSC, while Keimer et 



on YBCO compound, respectively. Their result of oblique lattice is not consistent with our 
result of triangular lattice. Both the single vortex structure and the vortex lattice structure is 
also different from theoretical predictions based on the dx2_y2-s admixture scenario. IBi'0'00 
The single vortex structure can have a significant influence on the structure of the vortex 
lattice in many- vortex problem. It is quite interesting to note that the four-fold symmetry 
of d-s admixed vortex is an implication of the mixed gradient coupling in its model. In our 
model, there is no mixed gradient term up to quadratic order. Very recently, Ichioka et al^ 
argued that a non-local correction is required and associated GL theory should enclose the 
fourth order mixed gradient term. They observed that the dx2-y2-wa:Ve component has the 
four-lobe shape and d^y-wave component the shape of octofoil. Although in this present 
paper we put emphasis on the effect of mixing through the direct coupling to the magnetic 
field, the result of Ichioka et al. strongly suggests that the dx2-y2-dxy admixture can leads to 
an oblique vortex lattice. Furthermore, even in the present GL model, for the intermediate 
field regime off the extreme regions near to Hd or Hc2, the order parameter distribution is 



V. DISCUSSIONS 



a/.0 and Maggio-Aprile et a/.0 observed oblique latt 



ices in their SANS and STM experiments 
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not spherically symmetric. Therefore further study is required for that region, because the 
non-spherical symmetry of the supercurrent distribution around a vortex is not compatible 
with the triangular lattice symmetry. 

The strong temperature dependence of Kes{T) in Eq. ( pT] ) reminiscent of i^2{T) in the 
microscopic consideration of conventional superconductors by Makiii and Eilenberge r§. 
Roughly speaking, the difference between H2{T) and k, for that case was due to the non- 
locality of the electromagnetic response of the superconductors. In our case, the strong 
T-dependence of K,es{T) has an interesting interpretation. In contrast to the low field limit, 
?/'_ is order of tp+ at high fields {H ~ Hc2) and low temperatures for z/p ~ 1, and the as- 
sociated paramagnetic current significantly affects the magnetization curve 4nM{H) in low 
temperature region. Several experiments reported the strong temperature dependence of 
the slope AndM/dH near H^2iT) in HTSC.§ 

VI. CONCLUSIONS 

We investigated a GL theory for vortex structures and magnetization in a d-wave super- 
conductor. In the GL theory, we assumed the d^^-yi-wave symmetry of the superconducting 
ground state, and the admixture of (i^2_j^2-wave and d^y symmetry in the presence of the vor- 
tices. The structure of an isolated single vortex was studied asymptotically and numerically 
at low field region {H ~ Hd). The isolated single vortex is similar to conventional s-wave 
vortex, and has almost spherically symmetric supercurrent distribution around it. The vor- 
tex lattice structure and magnetization were studied analytically at high fields near the 
upper critical field Hc2- The vortices arrange in a triangular lattice, and the magnetization 
curve 47rM(if) shows a strong temperature dependence for z/p ~ 1 due to the paramagnetic 
current effect. Some physical implications of the results were discussed. The results were 
also compared with the experimental observations and with those of d-s scenario. It was 
recognized that further study in the intermediate field region is valuable. 
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TABLES 

TABLE I. Characteristic lengths, scales of physical quantities, and phenomenological parame- 
ters in the GL theory in question. 



Characteristic lengths 


4m+ a+ ' 4m_Q;_' ''^ Sttc^ ^'oP' 


(In tlT*5» r'+flTl G^l P Tl ^il H ^ 




Fundamental parameters 


*2 = |a+|//3+, $„/27r = hc/le, k = \+/^+ 


Auxiliary parameters 


S = , T = m-3-/m4-3+, u = m-/m+, 

/ -J- 9 

X+ = Px/P, X- = MX+ 

C+ = M/3, C- = K+ 


Reduced units 


*±/*o = V'i, r/Te ^ T, r/^+ ^ r, 

B/F^ ^ B. A/(__f^ ^ A. J/(ci?^/4<_) ^ J 


^ Here H°2 is not the true upper critical field. 


The upper critical field can be substantially enhanced 



in this model. 
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FIGURES 

FIG. 1. Distribution of the magnitudes of the order parameter and the magnetic induction 

associated with an isolated single vortex, for different values of T and Vj^: T = 0.1,0.9 and 
u+ = 0.1,0.9. Other parameters are set as /x = 1.0, S(0) = 1.0, T = 1.0, x+ = 

FIG. 2. Temperature dependence of the k^q for three different values of i/p. Other parameters 
are set as /x = 1.0, S(0) = 1.0, T = 1.0, x+ = 0-8. 
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